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Abstract 

A brief survey of the theory of soliton perturbations is presented. The focus is 
on the usefulness of the so-called Generalised Fourier Transform (GFT). This is a 
method that involves expansions over the complete basis of "squared solutions" of 
the spectral problem, associated to the soliton equation. The Inverse Scattering 
CNJ ■ Transform for the corresponding hierarchy of soliton equations can be viewed as 

a GFT where the expansions of the solutions have generalised Fourier coefficients 
given by the scattering data. 
' The GFT provides a natural setting for the analysis of small perturbations to 

an integrable equation: starting from a purely soliton solution one can 'modify' the 
soliton parameters such as to incorporate the changes caused by the perturbation. 

As illustrative examples the perturbed equations of the KdV hierarchy, in partic- 
ular the Ostrovsky equation, followed by the perturbation theory for the Camassa- 
Holm hierarchy are presented. 
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1 Introduction 



Integrable equations are widely used as model equations in various problems. The inte- 
grability concept originates from the fact that these equations are in some sense exactly 
solvable, e.g. by the inverse scattering method (ISM), and exhibit global regular solu- 
tions. This feature is very important for applications, where in general analytical results 
(first integrals, particular solutions) are preferable to numerical computations, which are 
not only long and costly, but also intrinsically subject to numerical error. In a hydrody- 
namic context, even though water waves are expected to be unstable in general, they do 
exhibit certain stability properties in physical regimes where integrable model equations 
are accurate approximations for the evolution of the free surface water wave cf. [TJ I2T] . 

There are situations however where the model equation is not integrable, but is some- 
how close to an integrable equation, i.e. can be considered as a perturbation of an inte- 
grable equation. In such case it is still possible to obtain approximate analytical solutions. 
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There are several approaches treating the perturbations of integrable equations. One pos- 
sibility is to consider expanding the solutions of the perturbed nonlinear equation around 
the corresponding unperturbed solution and to determine the corrections due to pertur- 
bations. In other words, one represents the solutions u(x,t) in the form: 

u(x, t) = u(x, t) + Au(x, t), 

where u(x,t) is the solution of the corresponding unperturbed nonlinear evolutionary 
equation and Au(x,t) is a perturbation. The strength of the perturbation is measured 
by a parameter e, Au(x,t) = 0(e). By small (weak) perturbation one means 0<t< 1. 
Such perturbations can be studied directly in the configuration (coordinate) space, while 
the effect of the perturbations on the corresponding scattering data can be studied in 
the spectral space (usually the complex plane of the spectral parameter) of the associated 
spectral problem. 

For a direct study of soliton perturbations, one can use the multi-scale expansion 
method [29, 30j, introducing multiple scales, i.e. transforming the independent time 
variable t into several variables t n , (n = 0, 1, 2, . . . ) by 

t n = e% n = 0,1,2,..., 

where each t n is an order of e smaller than the previous time t n _i- Then, the time-derivative 
are replaced by the expansion (the so-called "derivative expansion") with respect to the 
multiple scales: 

oo 
n=0 

The dependent variable is expanded in an asymptotic series 

oo 

u(x,t) = ^2e n u n (x,t). 

n=0 

These expressions are substituted back into the equation, giving a sequence of equations 
for u n (x, t), corresponding to each order of e (each time scale t n = et). Solving the system 
of equations for u n (x,t), one has to ensure that there are no singularities in the solutions 
(i.e. that the solutions do not blow up in time, etc.). This may lead to some additional 
conditions on the functions u n (x,t) (or on the parameters in them), known as secular 
conditions. 

Several authors had used various versions of the direct approach in the study of soliton 
perturbations: D. J. Kaup JHS] had used a similar approach for the perturbed sine-Gordon 
equation. Keener and McLaughlin |69j had proposed a direct approach by obtaining the 
appropriate Green functions for the nonlinear Schrodinger and sine-Gordon equations. 
For a comprehensive review of the direct perturbation theory see e.g. [29], HI] and the 
references therein. 

In the spectral space, the study of the soliton perturbations is based on the perturba- 
tions of the scattering data, associated to the spectral problem. Such methods are used by 
a number of authors, for studying perturbations of various nonlinear evolutionary equa- 
tions, like the sin-Gordon equation [65], the nonlinear Schrodinger equation (TO], [38], 137] 
and of course, the KdV equation which is discussed in details in the following sections. 
The method is based on expanding the 'potential' (i.e. the dependent variable) u(x, t) of 
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the associated spectral problem over the complete set of "squared solutions", which are 
eigenf unctions of the corresponding recursion operator. 

The squared eigenfunctions of the spectral problem associated to an integrable equa- 
tion represent a complete basis of functions, which helps to describe the Inverse Scattering 
Transform for the corresponding hierarchy as a Generalised Fourier transform (GFT). The 
Fourier modes for the GFT are the Scattering data. Thus all the fundamental properties 
of an integrable equation such as the integrals of motion, the description of the equations 
of the whole hierarchy and their Hamiltonian structures can be naturally expressed mak- 
ing use of the completeness relation for the squared eigenfunctions and the properties of 
the recursion operator. 

The GFT also provides a natural setting for the analysis of small perturbations to an 
integrable equation. The leading idea is that starting from a purely soliton solution of a 
certain integrable equation one can 'modify' the soliton parameters such as to incorporate 
the changes caused by the perturbation. There is a contribution to the equations for the 
scattering data that comes from the GFT-expansion of the perturbation. 

In this review article we illustrate these ideas with several examples. Firstly we con- 
sider the equations of the KdV hierarchy and the KdV perturbed version - the Ostrovsky 
equation. Then we present the perturbation theory for the Camassa-Holm hierarchy. 

2 Basic facts for the inverse scattering method for 
the KdV hierarchy 

2.1 Direct scattering transform and scattering data 

The spectral problem for the equations of the KdV hierarchy is [82l 08] 

-y xx + u(x)y = k 2 m, (i) 

in which the real- valued potential u(x) is taken for simplicity to be a function of Schwartz- 
type: u(x) G <S(M), k G C is spectral parameter. The continuous spectrum under these 
conditions corresponds to real k. The discrete spectrum consists of finitely many points 
k n = iK n , n=l,...,N where real. 

The Jost solutions for (CQ) are as follows: f + (x,k) and f + (x,k) are fixed by their 
asymptotic when x — > oo for all real k ^ |82j : 

lim e~ ikx f + (x,k) = 1, (2) 

X— >QO 

f~{x, k) and/ - (re, k) fixed by their asymptotic when x — > — oo for all real k ^ 0: 

lim e ikx f~{x,k) = 1, (3) 

X— > — oc 

Since u(x) is real then 

f + (x,k) = f + (x,-k), and f-(x,k) = f-{x,-k). (4) 
In particular, for real k ^ we have: 

f ± (x,k)=f ± (x,-k), (5) 
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and the vectors of the two bases are related @: 



f~(x, k) = a(k)f + (x, -k) + b(k)f + (x, k), Im k = 0. (6) 

The coefficient a(k) allows analytic extension in the upper half of the complex fc-plane 
[H2] and 

a{k) = a{-k), b(k) = b(-k), \a{k)\ 2 - \b{k)\ 2 = 1. (7) 

The quantities TZ ± {k) = b(±k)/a(k) are known as reflection coefficients (to the right with 
superscript (+) and to the left with superscript (— ) respectively). It is sufficient to know 
TZ ± (k) only on the half line k > 0, since from ([7]): 7?. ± (— k) = 1Z ± (k) and also (171) 



\a(k)\ 2 = (1 - 1^(^)1 
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Furthermore TZ ± (k) uniquely determines a(k) [82J. At the points K n of the discrete spec- 
trum, a(k) has simple zeroes i.e.: 

a(k) = (k — in n )a n + -{k — m n ) 2 d n H , (9) 

The dot stands for a derivative with respect to k and d n = d(m n ), d n = a(tK n ), etc. The 
following dispersion relation holds: 

a{k) = exp ( i r mw) n ^ 

At the points of the discrete spectrum /~ and / + are linearly dependent: 

f~(x, in n ) = b n f + (x, iK n ). (11) 

In other words, the discrete spectrum is simple, there is only one (real) linearly indepen- 
dent eigenfunction, corresponding to each eigenvalue in n , say 

fn( x ) = /~0M«n) (12) 

From (fT2l and (J2J), (jSJ) it follows that f^{ x ) falls °ff exponentially for x — > ±oo, which 
allows one to show that f n (x) is square integrable. Moreover, for compactly supported 
potentials u(x) (cf. ( TTTT) and (jBJ)) 

b n = b(iK n ), b(-iK n ) = -—. (13) 

b„ 

One can argue [82J, that the results from this case can be extended to Schwarz-class 
potentials by an appropriate limiting procedure. 
The asymptotic of /~, according to (jSJ), (T5]), (TTTT) is 

f~(x) = e KnX + o(e KnX ), x -* -oo; (14) 
f~(x) = b n e- K " x + o(e- K " x ), x-^oo. (15) 



'According to the notations used in [S2] f + (x, k) = "4>{x, k), f (x, k) = <p(x, k). 
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The sign of b n obviously depends on the number of the zeroes of /~. Suppose that 
< Ki < k>2 < ■ ■ ■ < k>n- Then from the oscillation theorem for the Sturm-Liouville 
problem [3], /~ has exactly n — 1 zeroes. Therefore 

6n = (-l) n - 1 |6n|. (16) 

The sets 

S ± = {TZ ± (k) (k>0), n n , B$ = ^, n=l,...N} (17) 

ia n 

are called scattering data. Clearly, due to ffTUj) each set - S + or S~ of scattering data 
uniquely determines the other one and also the potential u(x) 

2.2 Generalised Fourier Transform 

The recursion operator for the KdV hiererchy is 



l„o , . 1 



4 y ' 2 



X 



L± = --d 2 + u(x) - - / dxu'(x) ■ . (18) 



±oo 



The eigenfunctions of the recursion operator are the squared eigenf unctions of the spectral 
problem: 

F±(x, k) = k)f, Ff{x) = F{x, iK n ), (19) 

where F^(x) are related to the discrete spectrum k n = in n . Using (fj]) and the asymptotics 
©, © one can show that 

L ± F± (x, k) — k 2 F ± (x, k) L±F± (x) = k 2 n F± (x) . (20) 

The Generalised Fourier expansion can be formulated as follows: 

Theorem 2.1 Suppose that f + and f~ are not linearly dependent at x = 0. For each 
function g{x) G <S(M) the following expansion formulas hold: 

1 f°° N / 

!!(■'■) = ±2rij jHk)F±(x,k)dk T J2 (atjFU x ) +9&F±{x] 

where Fj (x) = [-§^F ± (x,k))k=k j an d the Fourier coefficients are 

1 f°° 
~9 ± {k) = 7^y(^ FT )' where (9,F) = J g{x)F{x)d 

9 ^ = fc^f^'^' 



± 



9% 



kjdj 
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Proof: The details of the derivation can be found e.g. in [33), SH] - B 
In particular one can expand the potential u(x), the coefficients are given through the 
scattering data 



2 f°° ^ 
u(x) = ±— kn ± (k)F ± (x,k)dk + J2<iik j RfF±(x). (21) 

The variation Su(x) under the assumption that the number of the discrete eigenvalues is 
conserved is 

1 f°° N r i 

5u(x) = — 6KHk)F±(x,k)dk±2j2\Rpk j F±(x)+6RfF± . (22) 

71 J -°° 3=1 

An important subclass of variations are due to the time-evolution of u, i.e. effectively 
we consider a one-parametric family of spectral problems, allowing a dependence of an 
additional parameter t (time). Then 5u(x,t) = u t 5t + Q((5t) 2 ), etc. The equations of the 
KdV hierarchy can be written as 

u t + d x Q{L±)u{x,t) = 0, (23) 

where Q(k 2 ) is a rational function specifying the dispersion law of the equation. The 
substitution of (l22l) and ( !2T|) in ( |23l) due to ( !20l) gives a system of trivial linear ordinary 
differential equations for the scattering data: 

TZf ± 2tkn(k 2 )TZ ± = 0, (24) 

Rf t ±2tk j n(k 2 j )Rf = 0, (25) 

k ht = 0. (26) 

The KdV equation 

u t - Quu x + u xxx = (27) 

can be obtained for Q(k 2 ) = —4k 2 . 

Once the scattering data are determined from f[2~4"j) - fl2"B"l) one can recover the solution 
from (|2"2"]) . Thus the Inverse Scattering Transform can be viewed as a GFT. 



3 Perturbations of the equations of the KdV hierar- 
chy 

Let us now consider a perturbed equation from the KdV hierarchy, i.e. an equation of 
the form 

u t + d x n(L±)u(x, t) = P[u], (28) 

where P[u] is a small perturbation, which is also assumed in the Schwartz-class. The 
perturbed equation is, in general, non-integrable. One can expand P[u] according to 
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Theorem 12.11 and ut and u according to ( 1221) and ( 12T1) . Their substitution in ( 1281) now 
apparently leads to a modification of (124)) - (126]) as follows: 



Kf ± 2ikVt{k 2 )TZ ± 



R% ± 2 /7,; ; <)(/,•;)/.'; 



~2ka 2 (k) 
1 

'2kjoj 
2k j a 2 j Rf ' 



1 a,- 



+ 



(29) 
(30) 

(31) 



Note that due to ( 1311) as a result of the perturbation the discrete eigenvalues are time- 
dependent. Another feature is the contribution from the continuous spectrum: even if 
one starts with a pure soliton solution of the unperturbed equation (JZ ± (k,0) = 0) then, 
in general ^(k, t) ^ due to (|2H]). 

For practical applications it is easier to work with an equation for b n instead of ( 1301) . 
Such an equation can be obtained as follows. We notice that 



R 



ii.t 



R 



n,t 



1 

ib n h n j i 



—b n}t + b n ( — 
ia„ \ ia 



1 



-b n ,t 



n/ t 
1 ( 1 



%Qj n b n b Ti 



n/ t 



thus b n>t 
have 



^{Rt, t - b 2 n R~ t ). Then using (JUD and the fact that F~ = b 2 n F+, cf. (JTTJ we 



b n>t + 2ik n Q(k 2 l )b n 



P 6 2 F H 



Ak„L 



(32) 



As an example let us consider the adiabatic perturbation of the one-soliton solution of 
the KdV equation. The one-soliton solution is 



u s (x, t) = — 2/t 2 sech 2 ,2, z = Ki(x — £), £ = An\t + £ - 



The eigenfunctions are 



(33) 



The perturbed solution is 

u(x, t) = — 2K 2 sech 2 z + u r (x, t), z = Ki(t)[x — £(£)]. 

Here u r (x,t) is the contribution from the continuous spectrum (radiation). 
From (1311) we have 



Kit 



P[u s (z)]sech. 2 zdz . 



(34) 



(35) 



(36) 
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Writing bi = e 2 " 1 ^® and using ( 1361) and ( |32i) we obtain 

£ t = Ak\ - -^r / ( z + \ sinh 2z ) sech^dz. (37) 



For the reflection coefficient (I2"9"j) gives 



nt-8tk 3 7Z + = — / P[ Ms (2)]e- 2ife/K1 (A;-m 1 tanh2) 2 dz. (3* 



then according to [59] using approximations in Gelfand-Levitan-Marchenko equation one 
can obtain 
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Ur{Xj t ) = *± r n + {k)e ^ + 2ik Z / K1 f k + i Kl tanhz \ dL (3Q) 



71 dz \ k + IKi 

Alternatively, from (12111 it follows that 

u r (x,t) = — / kTZ + (k)F + (x 1 k)dk 

= 1 [°° kU + (k)e 2 ^ +2 ^ ( fc+ , mitanh ^ 2 dk. (40) 

Both formulae give an approximation of the radiation component since the z- derivative 
of tanh z can be neglected [76] . 

The perturbation results for the Zakharov-Shabat (ZS) type spectral problems have 
been obtained firstly in [H3] arid for KdV in |59j . As it has been explained, the perturbation 
theory is based on the completeness relations for the squared eigenf unctions. For the 
Sturm-Liouville spectral problem such relations apparently have been studied as early as 
in 1946 [1] and then by other authors, e.g. [7T1 148] . The completeness relation for the 
eigenfunctions of the ZS spectral problem is derived in [64J and generalisations are studied 
further in [Ml EH [Ml SD EZj , see also [ID] . 

Example: Ostrovsky equation. This equation has the form 



u t + u xxx - 6uu x = 7<9 V (41) 

where <9 _1 is an operator such that (d~ l u) x = u, in general not uniquely determined. 
The Ostrovsky equation can be viewed as a time-dependent nonlocal perturbation of the 
KdV equation fl2"Tj) . Here 7 is a constant parameter. The equation is often called the 
Rot at ion- Modified Korteweg-de Vries equation. It describes gravity waves propagating 
down a channel under the influence of Coriolis force. In essence, u in the equation can be 
regarded as the fluid velocity in the x— direction. The physical parameter 7 measures the 
effect of the Earth's rotation. More details about the Ostrovsky equation can be found 

e.g. in [Sam eg Eg. 

In the perturbation theory 7 < 1 plays the role of a small parameter. In order to 
ensure that the perturbation is decaying fast enough at x — ±00 we take the one-soliton 
KdV solution in the form 

u s = 2/t 2 /sinh 2 z (42) 
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which can be obtained formally from ( 1331) for ki^q = m/2. It is not continuous at z = 
but decays fast enough at x = ±00. Using the fact that 

^ = -sk + M « = + - 9( - 2) i (43) 

we obtain 

P[ Us ) = 1 d~ 1 u s = -27/ci [coth z - 0(z) + 0(-z)] , (44) 



which is an odd function of z and then (13611 gives K\,t — 0. Thus the amplitude of 
the soliton does not change under this perturbation. The computation of ( 1371) gives a 
correction to the velocity of the soliton: 

4 Conservation laws and perturbed soliton equations 



It is well known [881 E2] that the KdV equation is a completely integrable Hamiltonian 
system and possesses infinitely-many integrals of motion. These integrals can be con- 
structed from the scattering coefficients a(k) of the associated spectral problem (PQ) and 
are polynomials of the dependent variable u(x, t) and its x- derivatives: 

Pn{ 

U, U x , U X xy ■ ■ ■ 

) dx, (45) 

where P n is a polynomial with respect to u and its derivatives. Since a(k) is time- 
independent, it can be viewed as generating functional of integrals of motion |82j : 

00 T 

= < 46 > 

Skipping the details (see e.g. [8"2]). we provide here only the list of the first few integrals 
of motion: 

1 f°° 

h = — / u(x)dx; (47) 



h = ~ I u(x) 2 dx; (4E 



-00 
00 



-00 
00 



h = ~\ I (u 2 x (x) + 2u\x)) dx- (49) 



—00 
00 



1 f°° 

h = ~- (u 2 xx - 5u 2 u xx + 5m 4 ) dx; (50) 
^ J -00 



The KdV equation (|27p can be written as a Hamiltonian system 



d SH 

dx 5u(x) ' 
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where the symbol S/8u denotes variational derivative. Moreover, ( |5"T1) can be further 
represented in its Hamiltonian form with a Hamiltonian H: 



u t = {u,H}. (52) 
where H = I 3 (l4*9j) . The Poisson bracket is defined as 

The first three integrals of motion (T4Tj) -( l49l) have the same interpretation for all members 
of KdV hirarchy: The first one, I± is related to the algebraic structure of the Poisson 
bracket (15111) : it follows from the presence of the operator d/dx in the Poisson brackets. 
The integral (T4*H]) has a meaning of a momentum. It is related to the translation invariance 
of the Hamiltonian. Since H [u(x + e)} — H[u(x)] = 0, the expansion of J (H[u(x + e)} — 
H[u(x)])dx in e about e = gives (note that u(x) = SP/Su) 

0= f^L^dx= [ 6H d 6P dx = {p,m = p t , 

J 5u(x) dx J 5u(x) dx Su(x) ' 



Consider now the perturbed KdV equation 

We will seek the integrals of motion for the perturbed equation Ik in the form Ik = 
Ik + AJfc, k — 1,2, ... . Here Alk can be viewed as a correction to the integrals of motion 
of the unperturbed equation (1271) coming from the perturbations P[u]. After a direct 
integration of (1251) . one gets: 



Ali = / P[u] dx. (54) 

J — oo 

Then, multiplying (12"7|) by u(x, t), and integrating leads to: 

/oo 
uP[u] dx, (55) 
-oo 

and so on. 

As an illustrative example we will take again the Ostrovsky equation ( HTl) . Due to the 
concrete choice of the perturbation in the right hand side of ( HIT) , the integrals in ( 1541) 
and (j55l) vanish, so the perturbations do not contribute to these integrals: the first two 
integrals of motion are the same as for the KdV equation. The nontrivial contributions 
of perturbations to the integrals of motion in the Hamiltonian J3 are: 

AI 3 = 7: [ (d^ufdx. (56) 
* J -00 

Note also, that there is no second Hamiltonian formulation for the Ostrovsky equation, 
compatible with the one given above, i.e. the equation is not bi-Hamiltonian - indeed 
( |4T1) is not completely integrable for 7 7^ 0, [7j. 
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5 Perturbations to the equations of the Camassa- 
Holm hierarchy 

Closely related to the KdV hierarchy is the hierarchy of the Camassa-Holm (CH) equation 
[6]. This equation has the form 

Ut — u xxt + 2uju x + 3uu x — 2u x u xx — uu xxx = 0, (57) 

where uj is a real constant. It is integrable with a Lax pair [6] 

*xx = (~+A(m + u;))# (58) 
% = + + (59) 

where m = u — and 7 is an arbitrary constant. 

Both CH and KdV equations appeared initially as models of the propagation of two- 
dimensional shallow water waves over a flat bottom. More about the physical relevance 
of the CH equation can be found e.g. in El EE1 EU E21 ED EH HHJ ED] . The paper [75] 
suggests that KdV and CH might be relevant to the modelling of tsunami waves (see also 
the discussion in [IS]). 

While all smooth data yield solutions of the KdV equation existing for all times, certain 
smooth initial data for CH lead to global solutions and others to breaking waves: the 
solution remains bounded but its slope becomes unbounded in finite time (see p~3|. 151 15]). 
The solitary waves of KdV are smooth solitons, while the solitary waves of CH, which are 
also solitons, are smooth if uj > [HI ES] and peaked (called "peakons" and representing 
weak solutions) if uj = [HI HH El E31 EZ] • Both solitary wave forms for CH are stable 

It could be pointed out that the peakons appear also as travelling wave solutions of 
greatest height (for the governing equations for water waves) , cf . [HI [12], [86] . 

In geometric context, the CH equation arises as a geodesic equation on the diffeo- 
morphism group (if uj — 0) [SI [121 EH El] and on the Bott-Virasoro group (if uj > 0) 

CH equation also allows for solutions with compactly supported m(x, t), [TOj, however 
u(x, t) looses instantly its compact support, whether uj 7^ [12] or u = [78] . 

The problem of perturbation of the CH equation arises when one deals with model 
equations that are in general non-integrable but close to the CH equation. A perturbation 
could appear for example when one takes into account the viscosity effect [HI]. Another 
possible scenario comes from the so-called '6-equation' [28l Wf\ 

m t + bu)u x + bmu x + m x u = 0. 

The 6-equation generalizes the CH equation and is integrable only for 6 = 2 (when it 
coincides with the CH) and 6 = 3 (then known as Degasperis-Procesi equation) [801 Wf\ 152] . 
Qualitatively, the DP equation exhibits most of the features of the CH equation, e.g. the 
infinite propagation speed for DP was established in [13]. In [7H] it is suggested that 
DP (as well as CH) might be relevant to the modelling of tsunami waves (see also the 
discussion in [18J). 
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The hydrodynamic relevance of the 6-equation is discussed e.g. in [56, 5lj. There- 
fore, the solutions of the 6-equation for values of b close to b = 2 can be analyzed in 
the framework of the CH-perturbation theory. We can represent the equation as a CH 
perturbation 

m t + 2uu x + 2mu x + m x u = (2 — b)(uju x + mu x ) = P[u], 
for a small parameter e = b — 2. 



5.1 Inverse scattering and generalised Fourier transform for the 
CH spectral problem 

The CH spectral problem ( 1581) can be handled in a way, similar to the one, already 
outlined. For simplicity we consider the case where m is a Schwartz class function, u > 
and m(x,0) + u> > 0. Let is introduce the notation q(x,t) = m(x,t) + to. Then one can 
show that q(x,t) > for all t [9]. Let k 2 = — -r — Xu, i.e. 

A(fc) = _!(*»+ 1). (60) 

The spectrum of the problem (1581) under these conditions is described in [9]. The con- 
tinuous spectrum in terms of k corresponds to k - real. The discrete spectrum consists 
of finitely many points k n = in n , n = 1, . . . , N where K n is real and < K n < 1/2. The 
continuous spectrum vanishes for u = 0, [22J. 

All results (J2D - (TTTI) remain formally the same with the exception of ( flOi) which now 
has the form [151 [16] 

.(*) = exp (-iak - -L f" M1-J^WI\,A (61) 



where 



oo 

N 



n ---- / (W^-l)dr 



^ Vl-2« n y tt7 4P + 1 



is one of the CH integrals of motion (Casimir). 

With the asymptotics of the Jost solutions and ( |58|) one can show that 

LiF^x, fe) = ^F±(x, k) L ± F±(x) = ^-F±(x), (62) 

where A n = X(iK n ); F are again the squares of the Jost solutions as in ( TT9l) and 

L± = (d 2 - l^Uqix) -2 f dym\y)\ (63) 

i±oo 

is the recursion operator. The inverse of this operator is also well defined. 
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The completeness relation for the eigenfunctions of the recursion operator is [T5] 



U! 



y/q(x)q(y) 
1 



{x-y) 



2m 



F-(x,k)F+(y,k) 
ka 2 {k) 



dk 



N 



n=l 



F-(x)F+(y) + F-(x)F+(y)- 



(In 



(64) 



Therefore F ± , F^ and F^ can be considered as 'generalised' exponents. Like in the KdV 
case it is possible to expand m(x) and its variation over the aforementioned basis, or 
rather the quantities that are determined by m(x) and Sm(x) 
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q(x) 



'2m 



2 -^^FHx,k)dk + J2 2 -t^F n ± 



n=l 



X 



(65) 



/ S^qJi) dy = — SK+WF^^dk 
J±oo 2ttz J_ x \{k) 

N r 1 /?± i 

n=l 

The expansion coefficients as expected are given by the scattering data and their varia- 
tions. This makes evident the interpretation of the ISM as a generalized Fourier transform. 
Now it is straightforward to describe the hierarchy of Camassa-Holm equations. To every 
choice of the function fl(z), known also as the dispersion law we can put into corre- 
spondence the nonlinear evolution equation (NLEE) that belongs to the Camassa-Holm 
hierarchy: 

^jT(V8A + nfe)(^-i)=o. (67) 

An equivalent form of the equation is 

q t + 2qu x + q x u = 0, u = )fl{L ± ) - l) . (68) 

The choice Q(z) = z leads to u = u and thus to the CH equation (15 7p . Other choices of 
the dispersion law and the corresponding equations of the Camassa-Holm hierarchy are 
discussed in [161 E3] . 

By virtue of the expansions ( 1651) and (15511 the NLEE (1571) is equivalent to the following 
linear evolution equations for the scattering data: 

Kf{k) =f ikVl^X- 1 )^^) = 0, (69) 
i2j t ±« n fi(A^ 1 )i^ = 0, (70) 
«n,t = 0. (71) 

The time-evolution of the scattering data for the CH equation (1571) can be computed from 
the above formulae for Q(z) = z, see also [13, fT5] . 
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5.2 Perturbation theory for the CH hierarchy 

Let us start with a perturbed equation of the CH hierarchy of the form 

q t + 2qu x + q x u = P[u], u = -ft (At) (i/^ - 1 



(72) 



where again, P[u] is a small perturbation, by assumption in the Schwartz-class. It is 
useful to write (1721) in the form 



2 



(^) t ck/ + ft(L ± )L/--l 



1 r P(y) 



dy. 



(73) 



With the completeness relation (1641) one can deduce the gereralised Fourier expansion for 
expressions, like the one on the right-hand side of (1721 

Theorem 5.1 Assuming that f + and f~ are not linearly dependent at x = and g(x) £ 
<S(R), the following expansion formulas hold: 



g(y) 



1 

2ni 

N 



dy= ± — I g^QFffaQdk 



(74) 



where Fj(x) = [-^rF ± (x, k)] k=k . and the Fourier coefficients are 



ka 2 (k) 



kjOj L 



The substitution of the expansions (1741) for -P[w], (16"5"1) and (j6U|) into the perturbed 
equation (173"]) gives the following expressions for the modified scattering data: 



i\{P,F*) 
T 2ko?(k) ' 
A 3 -(P,i7) 
2%d]it* 
R%-Rf\ j>t ± kMI/X^R* 

A,- 



72^ =F«fcft(l/A)7e ± 



2%dj 



(75) 
(76) 

(77) 



From (1771) we obtain the following for the coefficient bf 

Aj 



6j,i + K j ft(l/A i )6 i 



, . .P,tiFf-F7 
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The 'perturbed' solution for the hierarchy in the adiabatic approximation can be recovered 
from the following expansion for u(x) with the 'modified' scattering data keeping the 
unperturbed 'generalised' exponents: 



1 



oo 



fcfi(l/A(Jfe)) 



N 




u(x) = ± 



2ni 



n ± (k)F ± (x, k)dk + J2 



— oo 



n=l 



This formula follows from the second part of (|68|) and (|65|) . Note that for the CH equation 
Qg?D u = u. 

6 Discussions 

We have presented a review of some aspects of the perturbation theory for integrable 
equations using as main examples the KdV and CH hierarchies. 

In our derivations we used completeness relations that are valid only given the as- 
sumption that the Jost solutions / + and /~ are linearly independent at x = 0. The 
case when this condition is not satisfied is quite exceptional, however this is exactly the 
case when one has purely soliton solution [33], HE] . Then one has to take into account a 
nontrivial contribution from the scattering data at k = [16] and some of the presented 
results require modification. E.g. (I3~TI) should be [58J 



since P(z) is an odd function and the additional term does not contribute. The meaning 
of the condition ( 1791) is that no shelf is formed behind the soliton [581 H6j. The presence 
of shelf for KdV equation is observed e.g. under the perturbation P[u] = eu [58l [76] . 

The corrections to the conservation laws due to perturbations have been used in studies 
of the effects of the disturbance on the initial soliton [5"S1 155] , or as a correctness check of 
results obtained otherwise [67] . 

The evaluation of the perturbation terms for the CH hierarchy could be technically 
difficult due to the complicated form of the CH multisoliton solutions [9} [15]. However 
the limit u — > leads to the relatively simple peakon solutions. Therefore, using the 
presented general formulae one should be able to access the perturbations of the peakon 
parameters. 
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P[u s (z)](zsech. 2 z + tanhz + tanh 2 z)dz. 



(78) 



In the presented example with the Ostrovsky equation 




(79) 
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